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Abstract. The neutral massless scalar quantum field $ in four-dimensional space- 
time is considered, which is subject to a simple bilinear self-interaction. Is is well-known 

2 r, 

from renormalization theory that adding a term of the form — to the Lagrangean 

has the formal effect of shifting the particle mass from the original zero value to m 
after resummation of all two-leg insertions in the Feynman graphs appearing in the 
perturbative expansion of the 5-matrix. However, this resummation is accompanied 
by some subtleties if done in a proper mathematical manner. Although the model 
seems to be almost trivial, is shows many interesting features which are useful for the 
understanding of the convergence behavior of perturbation theory in general. Some 
important facts in connection with the basic principles of quantum field theory and 
distribution theory are highlighted, and a remark is made on possible generalizations 
of the distribution spaces used in local quantum field theory. A short discussion how 
one can view the spontaneous breakdown of gauge symmetry in massive gauge theories 
within a massless framework is presented. 
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1. Introduction 



The traditional starting point of perturbative quantum field theory is a classical 
Lagrangean C which can be decomposed into a free (solvable) part Cq and an interacting 
part C int which describes the interaction. These objects get quantized and ^-matrix 
elements or Greens functions are constructed with the help of the Feynman rules. In 
this paper, we focus mainly on the extremely simple case of a free scalar massless field 
$ fulfilling the wave equation 

= d„d^(x) = 0, (1) 
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which can be decomposed into a negative and positive frequency parts according to the 
representation (k° = \k\, kx = k^x^ = k°x° — kx) 

/d 3 k r -> i 
-j= a{k)e~ ikx + a\k)e ikx (2) 



'2|Jfe| 

with the distributional commutation relations for the creation and annihilation 

'operators' pQ 

[a(k), a\k')} = 8&(k - k'), [a(k),a(k')\ = [a f (£), a\k')\ = 0, (3) 
where the full Lagrangean describing the dynamics of the field shall be given by 
C = C + Li n t with 

A) = \d^&§, C mt = - l -m 2 e. (4) 

Actually, we will avoid the simple way by solving the massive Klein-Gordon equation 
(□ + m 2 )$(x) = 0, which directly leads to the corresponding free massive scalar field. 
The aim of this paper is to highlight in which sense the massive theory comes out as 
a perturbative limit of the massless case. This transition, although seemingly trivial, 
shows many interesting features which can be observed also in 'more realistic' theories. 
E.g., the positive energy spectrum of the massless theory is continuous, whereas the 
energy spectrum of the new massive theory displays a gap between the vacuum and 
the lowest one-particle energy state. An example for an exactly solvable theory which 
displays such a behavior is the Schwinger model, i.e. massless quantum electrodynamics 
in 2+1 space-time dimensions [21 IB]- 

The free massless field $ is an operator valued distribution in the sense that it must 
be smoothly averaged over some space(-time) region according to the formal expression 

$(/) = / $(x)f(x)d 4 x (5) 



in order to yield an operator which is densely defined in the corresponding Fock space. 
E.g., one might assume that $(a;) is an operator applicable to the Fock vacuum |0). 
Then a simple calculation using the relations given above shows that $(x)|0) can not 
be associated with a vector in Fock space. Within the famous Wightman framework [1], 
/ in eq.(|5]) is an element of the Schwartz space iS(R 4 ) or rapidly decreasing functions. 
There are several good reasons for the use of the Schwartz space in the Wightman 
formalism. The Fourier transform is a linear isomorphism on 5(R 4 ). Correspondingly, 
the Fourier transform acts also as a linear isomorphism on the dual space of tempered 
distributions iS'(R 4 ). Furthermore, iS(R 4 ) contains the test functions with compact 
support D(R 4 ) C ^(R 4 ), which are used to express the locality and causality properties 
of field operators. E.g., the causal property of the scalar field can be characterized by 
test functions /, gG£>(R 4 ) via 

[$(/), $(<?)] = V/,</ with supp(/) ~ supp(^), (6) 

where supp(/) ~ supp(g) denotes the fact that every element in the support of / is 
space-like separated with respect to every element in the support of g. 
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2. Resummation of the 2-leg insertions 

The theory defined by the Lagrangean eq. (HI) is ultraviolet super-renormalizable and 
contains only diagrams with a very simple topology. Vacuum diagrams have a polygonial 
structure, and will play no relevant role in the forthcoming discussion, i.e. we argue that 
they lead only to a phase factor in the S- matrix and can be 'divided away'. The $ 2 -term 
which defines the interaction leads to 2-leg insertions in the free massless propagator. In 
renormalization theory, such insertions are used to switch from the 'bare' particle mass 
to the 'physical' mass, but in our case, we will start from the purely massless theory 
and investigate the effect of including an interaction mass term. 

In the literature, a perturbative expansion of the S'-Matrix according to 

S = 1 + - j j d 4 x 1 ...d 4 x n T{H in t(xi)'Hint(x2) ■ ••• • Tiint(x n )}, (7) 

n=l n ' J 

is widespread, where T is the time-ordering operator and T~Ci nt is the Hamiltonian 
interaction density, which is given by TCi nt = —C int in our case. It must be pointed out 
that the perturbation series eq. ([7]) is formal and it is difficult to make any statement 
about the convergence of this series for a general quantum field theory. Furthermore, 
two problems arise in the expansion given above. First, the time-ordered products 

T n (x 1 ,x 2 , ...,x n ) = {-i) n T {Hi nt {xi)Hi n t{x 2 ) ■ ... • U int {x n )} (8) 

are usually plagued by ultraviolet divergences when calculated in momentum space 
according to the Feynman rules. Still, these divergences can be removed by 
regularization, such that the operator-valued distributions T n can be viewed as well- 
defined, already regularized expressions [5j El Ej. In the present case, one has 
T\(x) = — ^Tj- : $(x) 2 := i : d n t :, such that ultraviolet divergences are absent in 
the diagrams considered in this paper. The colons denote normal ordering. Second, 
infrared divergences are also present in eq. ([7]). This is not astonishing, since the T n 's 
are operator-valued distributions, and therefore must be smeared out by test functions 
in 5(R 4n ). One may therefore introduce a test function g(x) GiS(R 4 ) which plays the role 
of an 'adiabatic switching' and provides a cutoff in the long-range part of the interaction, 
which can be considered as a natural infrared regulator [5], [8] . The correct expression 
for the infrared regularized S-matrix is given by 

oo oo ^ « 

s (9) = ^2S n {g) = / d 4 xi...d 4 x n T n (a;i, ...x n )g(x 1 )-...-g(x n ), (9) 

n=0 n=l n ' J 

and an appropriate adiabatic limit g — > 1 must be performed at the end of actual 
calculations in the right quantities (like cross sections) where this limit exists. This 
is not one of the standard strategies usually found in the literature, however, it is the 
natural one in view of the mathematical framework used in perturbative quantum field 
theory. Performing the adiabatic limit is also necessary to restore the full Lorentz 
invariance of the theory. 
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The Feynman propagator of the free massless scalar quantum field is given by 

d 4 k e~ ikx 



A F (x) = -i(0|T($(a;)$(0))|0) 
i 1 i „ 1 1 



(2tt) 4 k 2 + iO 



4n 2 x 2 -iQ Ati^x 2 4t/ (x ^ (10) 
where T is the time-ordering operator, P denotes principal value regularization and S 
is the one-dimensional Dirac distribution depending on x 2 = = (x ) 2 — (x 1 ) 2 — 

(x 2 ) 2 — (x 3 ) 2 = x\ — x 2 . 

Usually, the resummation of the particle propagator is performed by a formal 
calculation. One observes first that the lowest order non-trivial contribution to the 
S'-matrix is generated by the interaction term Ti{x) = — nnl ; $(x) 2 ; = % : £ int ; ; which 
enters the S'-matrix at first order in m 2 as 

Si(9) = ^ J d 4 *i:<&(zi) 2 :s(*i). (11) 
At second order, one has after Wick ordering of the field operators 

4 / — %Tfl^ \ ^ f 

s 2{g) = 2j — ) / d 4 xid 4 x 2 iA F (xi - x 2 ) : $(xi)$(x 2 ) : g(x 1 )g(x 2 ) 
+ other Wick contractions 

TTt f 

= -i>~k~ I d 4 x 1 d 4 x 2 A F (x 1 -x 2 ) :$(xi)$(x 2 ): g{x 1 )g(x 2 ) 

+ other Wick contractions, (12) 
and for S n (g) one obtains, taking the permutation symmetry of the T n into account, 

S n {g) = -i>^2~ / d 4 x 1 ...d 4 x n A F (x 1 - x 2 )...A F (x n -i - x n ) 

: $(a;i)$(x n ) : g(xi)...g(x n ) + other Wick contractions. (13) 

Assuming that one may perform a adiabatic limits g(xt) — > 1 in the right variables, it 
becomes clear that the particle-particle transition amplitude is described after a Fourier 
transform by the resummed Feynman propagator 

1 m 2 m 4 1 . 

+ 7T, ■ 7^. + 7T, ■ .,„ +-= ,. 9 „ 9 ■ , ■ ( 14 ) 



k 2 + i0 (k 2 + i0) 2 (k 2 + i0) 3 " /c 2 - m 2 + zO 
However, the situation is not as simple as it seems. First, the geometric sequence eq. 
(fl4"j) does not converge in the sense of tempered distributions in 5'(R 4 ). Second, the 
terms ~ (k 2 + z0)~ n for n > 2 are not uniquely defined in a distributional sense. These 
terms are too singular for small momenta k. This situation can be compared to the case 
of scalar fields in 1 + 1 space-time dimensions. There, even the zero mass limit of the 
propagator (k 2 — m 2 +i0) _1 does not exist, since 1/k 2 is not a tempered distribution in 
two dimensions [9] . There is no free massless scalar field theory in two dimensions which 
does not violate the Wightman axiom of positivity [TOl [TT] after infrared regularization. 
Third, the geometric series eq. (fl4l does not converge for arbitrary k. We will investigate 
these defects in the following section. 



Resummation of mass terms in perturbative massless quantum field theory 



5 



3. Solutions of the n-fold iterated wave equation 

It is obvious that the expressions ~ (k 2 + i0)~ n for n > 2 are infrared divergent for 
k = are therefore need to be regularized in a way consistent with the basic principles 
of local quantum field theory. If this were not possible, the theory would be infrared 
non-renormalizable. 

In order to accomplish this task, we consider first commutation relations for 
massless scalar fields which define the positive and negative frequency massless Jordan- 
Pauli distributions 

A ± (x) = i{$ T (x), $±(0)] = i(0|[$ T (;r), $ ± (0)]|0) , (15) 

which have the Fourier transforms 

A ± (A;) = J d 4 xA ± (x)e ikx = ±2wiQ{±k Q )5{k 2 ), (16) 

where is the Heaviside distribution. The fact that the commutator 

[$(x), $(0)] = -%A + (x) - iA~(x) =: -iA(x) (17) 

vanishes for spacelike arguments (with x 2 < 0) due to the requirement of microcausality, 
leads to the important property that the Jordan- Pauli distribution A has causal support, 
i.e. it vanishes outside the closed forward and backward light cone such that 

supp A(x) CTU? + , V ± = {x | x 2 > 0, ±x° > 0} (18) 

in the sense of distributions. Furthermore, the Jordan-Pauli distribution has a simple 
representation in configuration space 

A(x) = —sgn{x°)S{x 2 ), (19) 
2n 

and solves the wave equation DA(x) = with the Cauchy data A(0,x) = and 
(<9 A)(0,£) =6^(x). 

A further crucial observation is the fact that one can introduce the retarded 
propagator A ret (x) which coincides with A(x) on V + — {0}, i.e. A ret (ip) = A(ip) holds 
for all test functions in the Schwartz space y?GiS>(R 4 ) with support supp(t^) C R 4 — V . 
In configuration space, A ret (x) is obviously given by 

A ret (x) = —Q(x°)5(x 2 ), (20) 

2n 

and as a special case of the edge of the wedge theorem [1] it is known that the Fourier 
transform of the retarded distribution A ret (k) is the boundary value of an analytic 
function r(z), regular in T + := R 4 + iV + . It is given by 

A-\k) = -jyr^o = -P± + i^ g n(k°)5(k 2 ). (21) 

The analytic expression for the Feynman propagator is recovered from the observation 
that 

A F {x) = -A re \x) + A" (re), (22) 
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(0|T($(a;)$(0))|0) = e(x°)[$(x), $(0)] + (0|$(0)$(x)|0), (23) 
and consequently 

A F (k) = -A^(k) + A-(k) = (24) 

The retarded distribution A ret (x) is a weak solution of the inhomogeneous wave equation 
UA ret (x) = 5 {4) {x) = —DAp{x). (25) 

Note that the Feynman propagator can also be written as 

A F (x) = -e(x°)A + (x) + G(-x°)A-(x). (26) 

For later use we introduce now distributions E n (x), fulfilling the n-fold iterated 
wave equation D n E n (x) = 0, with the properties 

E x (x) = A(x), (27) 

DE n+1 (x) =E n (x), (28) 

E n {Xx) = X 2n ^E n {x), AGR, (29) 

satisfying the complete set of Cauchy data at x° = 

{d%E){0,x) = 0, k — 0, 2n - 1, (^ n " 1 f?)(0, x) = 5^{x). (30) 

The E n 's can be constructed in a straightforward way. We consider first the 
distribution 

E(x) = —sgn{x°)Q{x 2 ). (31) 
8n 

From 

d v Q(x 2 ) = 2x u Q\x 2 ) = 2xJ(x 2 ), (32) 

d u (2xj(x 2 )) = 85(x 2 ) + 4x 2 5'(x 2 ), (33) 

one obtains by means of the identity (x 2 5(x 2 ) = 0) 

x 2 5'(x 2 ) = -f-(x 2 5(x 2 )) - 5(x 2 ) = -5(x 2 ) (34) 
ax 2 

DQ(x 2 ) = 45(x 2 ) (35) 

and it is a simple exercise to show that also 

□ (sgn(a; o )0(x 2 )) = 4sgn(a; )5(a; 2 ) (36) 
holds. Therefore, one may define 

E 2 (x) := E(x) = -L sgn (x )Q(x 2 ), (37) 

07T 
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since DE(x) = A(x) and E 2 has the correct scaling behavior and fulfills also all other 
requirements given above. From OE 2 (x) = A(x), one obtains in momentum space 

p 2 E 2 {k) = -± S gn(k°)5(k 2 ) (38) 
and by means of the identity eq. fl34l) with x replaced by k one may write 

E 2 (k) = ^sgn(k°)5'(k 2 ). (39) 

All these calculations are formal to a certain extent, but correct. However, one should 
bear in mind that the meaning of the distribution 5'(k 2 ) is rather defined by E 2 given 
above than vice versa. Note that the positive frequency part 

l E 2 (ky = ^e(k°)5'(k 2 ) (40) 

is not unambiguously defined. Only derivatives of E in configuration space can be split 
unambiguously into positive and negative frequency parts, which emerge, e.g., in the 
photon propagator in the Landau gauge 

fiv _ k»k" 

K U andau (k) ~ % 2+ fc y° , (41) 
in the term ~ k^k u /k 4 . We give here simply the regularized result for E 2 . One has 

E t( x ) = TT^log ^ , (42) 



or 



16vr 2 "°V A 2 , 



E 2 {x) = +-^log\x 2 /\ 2 z \ + -Ugn(x°)©(x 2 ) 

lb7l Z 107T 

i 1 
E+(x) = log |x 2 /A^| + sgn(x°)0(x 2 ) 



16vr 2 01 ' 16vr 
where A-r, is a renormalization length scale, and combining E 2 and E 2 indeed gives 

E+(x) + E 2 (x) = ^sgn(x°)Q(x 2 ). (43) 

07T 

This result is directly related to the fact that the convolution of the Feynman propagator 
in configuration space appearing in 5*3 which corresponds to the formal expression 
~ (k 2 + i0)~ 2 in momentum space (see eq. [T3"j) leads to an integral of the form 

T^(k) = Tz\ f 1 

2lJ [J (2tt) 4 x\ -i0(x- xx) 2 - i0_ 

4 ( 2V og (-¥) ^Tik^ 1 ^-^^' (44) 

when the adiabatic limit in the variable X\ is performed. 1Z denotes the regularization 
procedure with renormalization scale \n- Obviously, the regularized expression for 
T^(x) is defined up to a constant in real space or up to a local distribution in momentum 
space ~ S^(k). The scaling symmetry of the Feynman propagator A^( Ax) = A~ 2 A^(x), 
A G R, is spontaneously broken by regularization, such that a corresponding scaling law 
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r^-(Ax) = T^(x) does not hold as one might expect naively from the formal definition 
of r 2 (x) in eq. f|T3|) . The adiabatic limit does not exist for the S n for n > 3 for the 
C-number distributions. 
Finally, from 

U{x 2 ) n = An{n + l){x 2 ) n -\ (45) 
we immediately deduce 

£ "W = (2 T )4 n - 1 ( n -l)!( n -2)! ^»"' 2s ^ » e ^»' < 46 » 
4. Yang-Feldman equations 

Since the E^s introduced in the previous section all have causal support, we can define 
retarded distributions 

E r f(x) = — r- 1 w -(x 2 ) n - 2 e(x )&(x 2 ), n>2, (47) 

n v ; (27r)4™- 1 (n- l)!(n-2)! v ; \ ) \ h - > v ; 

fulfilling □™£ , ^ et (x) = 8^> (x). Roughly speaking, the distribution E^ et corresponds to 
the expression (— l) n (k 2 + ik°0)~ n in momentum space (the minus signs stem from the 
fact that the wave operator □ corresponds to —k 2 in momentum space). Note that 
the E^ et still respect the scaling law eq. (1291) . Performing now a resummation of the 
distributions in configuration space corresponding to minus 

-i 2 4 

1 m m 

k 2 + tk°0 + {k 2 + tk°0) 2 + {k 2 + ik°0) 3 + "' ^ ' 
instead of minus 

1 m m , , 

+ ,,2 I 77^ + 7U I ,nvs +-. ( 49 ) 



leads to 



fc 2 + «0 (P + iO) 2 (P + iO) 3 

oo 

5])(-l) n - 1 (m 2 ) n - 1 S ; 



n=l 

,2\n-l 



(50) 



1 ( 2\r 

- A-W + E 4.->'")V-2)! '- 2 )"' 29 '-°) 9 '- 2 ' 

[ 2V X 2 V 

where J\ is the Bessel function of order 1. The result in eq. (1501) is indeed the massive 
retarded propagator A™*, i.e. the Fourier transform of —(k 2 — m 2 + z/c°0) _1 ! 

A way to introduce the resummed retarded distribution constructed above is 
supplied by the Yang-Feldman formalism. One may write the interacting scalar field as 
a formal power series $j n < = X^o( m2 ) n< ^™' ^ na ^ fulfills the equation of motion 

n<& int = -m 2 g§ inU (51) 
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where $o is the free massless (incoming) field. The higher terms are recursively defined 
by 

$„(x) = - / d 4 Xl A ret {x - x 1 )g{x 1 )$ n - 1 {x 1 ). (52) 



Without the test function g, the integral eq. (152]) would be meaningless even if one 
chooses a massive scalar field for $ . F° r example, the expression 



J d 4 xiA ret (x - xi)$(xi)g{x\) 



1 

47T 



d 3 xi— p — - , 1 " — — $(x u — \x — x\\, xi) (53) 



g(x° - 




X — X\\ 






X 


'- Xi 





becomes meaningless in the limit g — > 1 and contains a volume divergence, and the 
situation is similar in the massive case. Such issues are rarely discussed in the literature, 
and completely ignored in the original work of Yang and Feldman [12] . Some attempts 
to deal in a rigorous style with the problem of the adiabatic limit can be found in 

u tram. 

The inductive construction of the retarded distributions can also be understood 
from the following simple calculation in configuration space. Convolving A ret with an 
E^ et leads to an integral for n > 2 

(2tt)- 2 4 1 -" 
(n- l)!(n-2)! X 

d 4 xxe(x?)5(^)((x - Xl ) 2 ) n - 2 &(x - x?)6((x - Xl ) 2 ). (54) 



I x 



The integral / vanishes outside the closed forward light cone due to Lorentz invariance 
and the two 0-distributions in eq. (|54|) . Therefore we can go to a Lorentz frame where 
x = (x° > 0,0) and (x — Xi) 2 = (x ) 2 — 2x°\x*i\ due to the 5-distribution in eq. (154]) . 
and we obtain 

i( X °>o,o)= 1 ; — — x 

(n — l)\(n — 2)1 
" >,:i ((x - Xifr-^Or ) 2 - 2x o |fi|)0(^ - |^i|) 



2|xi 

(2 7 r)- 1 4 1 - n 
fn - lWn - 2) 



d|fi||fi|(((x ) 2 - 2x |fi|) i ) n - 2 e(x72 - \xx\) 



(2vr)- 1 4 1 -« x ((x ) 2 )"-^ 



(n-l)!(n-2)! 4n(n-l) 
and the Lorentz invariant expression for J becomes 

/(*) = (^4^ (x2)n " 0(X ° )0(x2) = (56) 
i.e. we recover E r ^ x as defined above. This shows that a formalism using retarded 
products of distributions has an advantage compared to the common strategy to work 
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with time-ordered products. It is clear that the adiabatic limit for the final integral 
eq. fl53l) defining the Yang-Feldman field operator cannot be performed, however, 
the adiabatic limit exists for the full recursively defined retarded propagator and this 
allows to reconstruct the massive theory. Of course, according to Haag's theorem [15 J, 
representations of the canonical commutation relations algebra to different masses are 
inequivalent, and an attempt to express the massive field operator by the massless field 
would be futile. 



5. Relation to the operator product expansion 

In this section, we briefly describe the relation of our findings to a well-known technique, 
namely the operator product expansion (OPE), which goes back to a work of Wilson 
[IB] . The OPE provides a method to expand singular operator products as a sum of 
nonsingular local operators with coefficients being (singular) C-number distributions. 
The problem with the series eq. (TBI) which is solved within the framework presented 
above can be related to the fact that only quantities defined at short distances compared 
to the characteristic length of the physical system under consideration, i.e. short 
compared to the length scale 1/m in our case, can be expanded in powers of m 2 . 

At short space-time distances, one may consider the time-ordered product of two 
scalar fields and expand the corresponding vacuum expectation value 

(0|T($(x)$(0))|0) = C 1 (x 2 ) + C m2 (x 2 )(0|[$ 2 ]|0) (57) 

with Wilson coefficients C\ and C m 2, where (0|[$ 2 ]|0) is the vacuum expectation value 
in the free massive theory of the regularized composite operator 

9 

Ttl 

(0|[$ 2 ]|0) = T ^log(m 2 /^) (58) 

with a renormalization scale fi-jz [IT] . Clearly, the whole propagator cannot be expanded 
in powers of m 2 due to the logarithmic dependence of (0|[$ 2 ]|0). However, the Wilson 
coefficients can be expanded in powers of m 2 , and the corresponding power series are 
even convergent. An explicit calculation, performed in Euclidean space as it is common 
in the literature, leads to (r 2 ~ — x 2 ) 

CAr 2 ) = — 5-^(1 + mrlAmr) log(yU^r/2) + even powers in mr), (59) 
4:7T z r z 

C m ,(r>) = 2J ^, (60) 
mr 

where I\ is the modified Bessel function of order 1. Note that the //^-dependence of 

d(r 2 ) is canceled by that of (0| [$ 2 ] |0). 

The correction to the commutator or retarded propagator originates from the cut 

of log(yU7^r) in the complex plane of r 2 . From eq. fl59l) one obtains 

I m / — 

7r — ——mrlximr) = —Ii(myr 2 ), (61) 

47HH 47rVr 2 
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and returning to Minkowski space by further using the identity I\(z) = \J\(iz) ) eq. (1BT1) 
is in perfect agreement with the corresponding Bessel term in eq. fl50l) . The problem 
of the series eq. ( 1T4"1) is related to the non-trivial m 2 -dependence of the local composite 
operator vacuum expectation value (0|[$ 2 ]|0) in the framework of the OPE. 

6. Mass terms for gauge fields 

It is an interesting question to what extent one could describe theories with interacting 
massive vector bosons on the basis of massless gauge fields. To make a step towards this 
direction, we consider a pure SU(N) Yang-Mills theory without fermions below, and call 
the massless gauge bosons gluons for short in the forthcoming. 
As a first order coupling, one can choose 

T x {x) = ig YM fabc{\ :A^ a {x)A vb {x)F^(x): - :A^{x)u b {x)d»u c {x):}, (62) 

where = d v A v a — d^A u a is the free field strength tensor, u a , u a are the (fermionic) ghost 
fields, and the f abc are the SU(N) structure constants corresponding to the color indices 
a,b and c. This first order coupling contains the well-known three-gluon vertex and a 
ghost-antighost-gluon coupling. A commutator term ~ gf abc A v h A^. ist absent in the free 
expression for F^, since this term contains an additional Yang-Mills coupling constant 
g YM , and the resulting four-gluon coupling term generated by such an additional non- 
abelian term would be of the order q 2 . 

In the 't Hooft-Feynman gauge, the asymptotic massless free fields in eq. fl62l) 
satisfy the commutation relations [H] 

iS ab g^A T (x - y) (63) 

-t5 ab A^(x-y), (64) 

and all other {anti-}commutators vanish. It can be shown that the coupling eq. ( 1621) 
respects perturbative gauge invariance which corresponds to the full BRST symmetry 
of the theory [T9J [20]. The introduction of ghost couplings is necessary to preserve 
perturbative quantum gauge invariance at first order and unitarity at second order in 
g YM . We give here a short definition of perturbative BRST symmetry based on a gauge 
charge operator Q (for details refer to [221 123])- 

Defining the gauge charge Q, which is the generator of gauge transformations, by 

Q : = J d 3 xd^ a (x) d u a {x), (65) 

where summation over repeated indices is understood, one is led to the following (anti-) 
commutators with the fields: 



[A$(x),A%\y)] = 

and 

{u^(x),u^(y)} = 



[Q, A£] = i&Ua, [Q, Ff] = 0, {Q, u a } = 0, {Q, u a } = -id,A». (66) 
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The operator Q has been introduced for the first time by Kugo and Ojima [21] • The 
pseudounitary gauge transformation of a field operator like, e.g., Aj, is implemented by 



A'^ = e~ iXQ A^e iXQ , AgR. (67) 
Calculating [Q,Ti(x)] leads after some algebra to the interesting result 

[Q,T x ]=id^ /v (68) 

T£ /x = igf abc {:u a A ub Fr: +~ :u a u b d»u c :}. (69) 



The first order coupling is gauge invariant through the presence of a ghost coupling 
term in T\ in the sense that the commutator of Q with T\ is a pure divergence in the 
the analytic sense, i.e. an infinitesimal gauge transformation adds only a divergence to 
the original interaction term. 

This actual definition of perturbative gauge invariance can be generalized to higher 
orders of perturbation theory by the requirement that [22] [23] 

n 

[Q, T n (xi, ...£„)] = i d^T^xt, ...x n ) = (sum of divergences), (70) 
i=i 

where T£„ is a mathematically rigorous (regularized) version of the time-ordered product 

T n ^(x 1 ,...,x n )' = 'T(T 1 (x 1 )...T 1 ') 1 (x / )...r 1 (x n )). (71) 

The identities eq. ( ITT]) imply the well-known Slavnov- Taylor identities, which express 
gauge symmetry on an analytic level for Greens functions [21]. 

As aforementioned, the four-gluon term ~ g^ M is missing in T\. This term appears 
in a natural way as a necessary local normalization term of the gluon-gluon scattering 
diagram at second order in order to preserve gauge invariance in the present framework. 
Q has the important property that it is nilpotent, i.e. Q 2 = 0. This basic property of 
Q and the so-called Krein structure on the Fock-Hilbert space [251 126] allows to prove 
unitarity of the S'-matrix on the physical Hilbert space Ti p h ysi which is a subspace of 
the Fock-Hilbert space T containing also the unphysical ghosts and unphysical degrees 
of freedom of the vector fields. 

The basic question is now whether one can add a mass term to the interaction, such 
that the gauge symmetry remains preserved in a certain 'hidden' way. Adding simply a 
Proca mass term to T\ 

T?' A = l -ml:A«AZ: (72) 

would destroy gauge invariance at first order, since 

[Qr.AlA£:] = 2i:dpV A£: (73) 

cannot be written as a divergence. But this defect can be remedied by adding 
additionally a ghost mass term 

! = im a :u a u a := -im a :u a u a : (74) 
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to Ti, since 



[Q, '■ u a u a :} =:{Q,u a }u a : 



:u a {Q,u a }:= i -.Uad^A^:, 



(75) 



and so 




vector boson and ghost mass is equal for a certain fixed color in the 't Hooft-Feynman 
gauge, and an analogous result holds for arbitrary renormalizable gauges. It is 
clear that a resummation of the mass term will change the theory completely, and 
massive longitudinal gluon modes will become physical states. This will necessitate the 
introduction of new physical and unphysical Higgs (ghost) fields and an adjustment 
of the boson masses in the theory according to the symmetry breaking scenario, in 
order to save the consistency of the new theory. The full story is therefore not so 
simple [27]. There original gauge symmetry of the massless theory leaves its traces in 
the massive theory, but one should not expect that is it possible to construct a Higgs 
free electroweak model of interacting massive vector bosons based on the considerations 
presented above. However, it should be pointed out that the widespread statement in the 
literature that the Higgs mechanism is responsible for the particle masses is misleading 
to some extent. One could also argue that the presence of Higgs fields is necessary under 
certain circumstances in order to save the consistency of the theory, when vector bosons 
become massive. 

7. Concluding remark 

The series eq. ( 1481) converges formally to the massive retarded propagator, but not in 
the distributional sense, i.e. not in the weak* topology. This can be realized in a simple 
manner by the observation that every individual term in eq. (HHl) displays a singular 
behavior on the light cone, whereas the massive propagator is singular on the mass 
shell. It is instructive to compare the situation to the following simple example in one 
dimension. We consider distributions 



where 5^ (x) denotes the k-th. derivative of the one- dimensional Dirac distribution here. 
Applying d n to a test function ip(x) G<S(R) results in 



and all d n have local support supp(d n ) = {0}. If the test function can be continued to a 
holomorphic function in an open region containing the real axis (i.e. ip is analytic in the 




(77) 





fc=0 
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sense of (p€.S(R) fl C U (R), and not only a rapidly decreasing smooth C°°(R) function), 
then the Taylor expansion eq. (1751) indeed converges to 

' p(*>(0) 
fc! 



lim = £ = ¥>(1) (79) 

fc=0 

and for the restricted space of test functions «S(R) fl C^(R) one obtains a weak limit 

fc! 



f— lr 

lim d„ = y^f^(z) = (y(z-l), supp(£(x-l)) = {l}. (80) 

fc=0 

This suggests the speculative question whether the choice of a modified space of test 
functions would have advantages compared to 5(R 4 ) concerning the consistency of 
the mathematical framework of non-perturbative quantum field theory. Some ideas 
to generalize the admissible test function spaces in quantum field theory have been 
investigated by Jaffe [28] ■ The idea to construct a 'hyperfunction quantum field theory' 
has been followed in [29l [30] . There, the main technical problem is the fact that rapidly 
decreasing holomorphic functions do not have compact support, such that one has to 
look for a new concept to express the fundamental causality condition which is expressed 
in the Wightman formalism by eq. ([6]). 
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